Quantum computers can search rapidly by using almost any selective transformations 
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The search problem is to find a state satisfying certain properties out of a given set. Grover's 
algorithm drives a quantum computer from a prepared initial state to the target state and solves the 
problem quadratically faster than a classical computer. The algorithm uses selective transformations 
to distinguish the initial state and target state from other states. It does not succeed unless the 
selective transformations are very close to phase-inversions. Here we show a way to go beyond 
this limitation. An important application lies in quantum error-correction, where the errors can 
cause the selective transformations to deviate from phase-inversions. The algorithms presented here 
are robust to errors as long as the errors are reproducible and reversible. This particular class of 
systematic errors arise often from imperfections in apparatus setup. Hence our algorithms offer a 
significant fiexibility in the physical implementation of quantum search. 

PACS numbers: 03.67.Ac, 03.67.Lx, 03.67.Pp 



I. INTRODUCTION 

Suppose we have a set of N items, j = 0, 1, 2, . . . A^ — 
1, and a binary function f(j) which is 1 if j satisfies 
certain properties (e.g. if it is a solution to a certain 
computational problem) and otherwise. Let T be the 
set of M items for which f{j) = 1, i.e. T = = 1} 

and \T\ = M. Consider the situation when the items 
are not sorted according to any property, but /(j) can 
be computed by querying an oracle that outputs f{j) for 
any input j . The search problem is to find an element of 
T (i.e. a solution) using the minimum number of oracle 
queries. The best classical algorithm for this problem 
is to randomly pick an item j, use an oracle query to 
check whether j G T, and then repeating the process till 
a solution is found. On the average, it takes 0{N/M) 
oracle queries to succeed, since M/N is the probability 
of the picked item to be a solution. 

In a quantum setting, Grover's search algorithm [Ij] 
provides a much faster way. The A^ items are encoded as 
basis states \ j) of an A^-dimensional Hilbert space, which 
can be realized using n = log2 N qubits (without loss of 
generality, we assume A" to be a power of 2). The ini- 
tial unbiased state is chosen as the equal superposition 
state, {\/\/N) generated by applying the Walsh- 

Hadamard transformation on |0). The target state \t) 
can be any normalised state X^jeT % within the target 
subspace, since measuring \t) will always give a solution. 
Grover's algorithm obtains \t) by applying 0[^J N /M) 
iterations of the operator Q — WloWIt on W\0). Here 
It = E,(-1)'^*|J)(J| and lo = are the 

selective phase-inversions of \t) and |0) states respec- 
tively. Grover's algorithm thus provides a quadratic 
speedup over the classical algorithm, as each iteration 
of Q uses one oracle query to implement It- 

Grover showed that his algorithm works even if the 



Walsh-Hadamard transform W is replaced by almost any 
unitary operator [/ In this case, the initial state C/|0) 
is a general (not necessarily equal) superposition of the 
basis states. The operator Q — UloWit is iteratively 
applied to U\0), and the target state \t) is obtained after 

0{l/au) iterations, where au = 



C/,oP is the 



projection of U\0) on the target subspace (for U — W, 
aw = y/M/N). As the probability of getting a target 
state upon measuring U\0) is afj, the target state can be 
obtained classically by 0(l/a^) preparations of J7|0) and 
subsequent projective measurements. Hence, the quan- 
tum algorithm provides a quadratic speedup over this 
simple scheme by doing the same job in 0{l/ai/) steps. 
This generalization is known as quantum amplitude am- 
plification 0, 01 , and forms the backbone of many other 
quantum algorithms. It has an important application 
when in a physical implementation W gets replaced by U 
due to some unavoidable error. The algorithm succeeds 
as long as U and W can be consistently implemented 
even when we do not know their precise form, making it 
intrinsically robust against certain types of errors. In the 
case of quantum search, provided ajj ^ aw, there is not 
much of a slowdown and hence almost any transforma- 
tion is good enough. 

Quantum amplitude amplification often fails, however, 
when the selective phase-inversions {It,Ia} are replaced 
by other selective transformations, say {StiSo}. Con- 
sider the simple case when St = Rf = J2j 
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and 5*0 = Rq = J2j ^'''^ are the selective phase 

rotations of \t) and |0) states by angles (j) and </? respec- 
tively. The well-known phase matching condition [j, Q 
demands |0 — (/sj ^ ajj for quantum amplitude ampli- 
fication to succeed. This is a very strict condition for 
au <^ 1, while the quadratic speedup is not of much use 
for large aij. In fact, systematic phase mismatching (i.e. 
|(/)— (^1 -j^C au) is known to be the dominant gate imperfec- 
tion in implementing quantum amplitude amplification, 
posing an intrinsic limitation to the size of database that 
can be searched l6il. 
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In this work, we show that a successful quantum search 
can be obtained with almost any selective transforma- 
tions {St, So}, provided their inverse transformations 
{5'/,S'q} are also available. This is useful in situations 
where the errors are reproducible (i.e. every time we 
ask for the transformation A the system implements the 
transformation B) as well as reversible (i.e. whenever we 
ask for the transformation the system implements the 
transformation B^). For instance, such systematic errors 
arise when there is incorrect calibration of the instru- 
mentation. In the following, we present two algorithms 
in this category, one iterative and the other recursive. 

In section II, we consider the case of diagonal selec- 
tive transformations, which rotate the phases of the de- 
sired states by any amount (unlike the selective phase- 
inversions that change the phase by tt) but leave all the 
other (non-desired) states unchanged. We then construct 
an operator which yields a successful quantum search al- 
gorithm when iterated on the initial state, and we show 
the algorithm to be optimal up to a constant factor. This 
iteratve algorithm does not work in the case when diago- 
nal selective transformations also perturb the non-desired 
states. In section III, we design a recursive quantum 
search algorithm for such transformations provided they 
are not too far off from the selective phase-inversions. 

The algorithm requires 0{l/a^'^^'^*'^°^) queries, where 
At = II 5*4-/(11 and Ao = HS'o — /oil are the distances of se- 
lective transformations from the corresponding selective 
phase-inversions, assumed to be small. It is straightfor- 
ward to extend the above two algorithms to situtations 
where the selective transformations are non-diagonal. We 
describe that in section IV, together with possible appli- 
cations of our algorithms to quantum error correction, 
quantum workspace errors and bounded-error quantum 
search. 



II. ITERATIVE ALGORITHM 

Consider those selective transformations {S't,S'o} 
which rotate the phases of the desired states by arbi- 
trary angles but leave all the other states unchanged. In 
case of |0), there is only one desired state and 5*0 = Rq — 



/-(I 



')|0)(0|. In case of \t), there can be multiple 



target states and the rotation phase can be different for 
different target states, so St = Rt = J2j ^ 
we iteratively apply the generalized quantum amplitude 
amplification operator Q = UR^U^ Rt on the initial state 
C/|0), we will not succeed in getting a target state unless 
the phase matching condition is satisfied. 

Instead, we iteratively apply a different operator, T — 
U Rq'^W R\U R'^W Rt, on the initial state C/|0). It uses 
two oracle queries, one for Rt and another for r\. It also 
uses /^Q^ = Rq"^ along with R^. Thus, unlike Q, it makes 
explicit use of the inverse transformations {/ij, Ob- 
serve that T is a product of two selective phase rota- 
tions: URq'^U^ is a rotation by —f of the state U\0), 



and R\u R^U'^ Rt is a rotation by Lp of the state R\U\0). 
We therefore have 



T ^UR^^U'^R'^, \a)=RlU\0). 



(1) 



Let |r) be a state orthogonal to \a) in the two- 
dimensional subspace spanned by U\0) and \a), such that 
up to an overall phase 

;7|0) = cos6i|cr) +sin6i|T) = cos6ii?jc/|0) + sin6l|r). (2) 

For a general vector |-!/)(a, 6)) = a\a) + 6|r) in this sub- 
space, we have 

T\tp{a,b)) = C/i?o'^[/^i?^|7/;(a,6)) = URg'^U^lijiae'f ,b)). 

(3) 

As UR^'^Wli^) = \t{j)-zU\0) withz^ {l~e-'^){0\U\'ilj), 
we have 

T\^{a,b)) = |?/'(ae*'^ - z cos 61, 6 - z sin 6*)). (4) 

Hence, T preserves this two-dimensional subspace. For 
any vector within this subspace, we can also write R'§ = 
e'^R-^, and so T is equivalent to UR^'^WR-'^ up to 
an overall phase, i.e. 



r = URa'^U^R-'P 



(5) 



The above operator is a special case of the general- 
ized quantum amplitude amplification operator with |t) 
as the effective target state. It satisfies the phase match- 
ing condition by construction. One may wonder that 
the phase matching condition is not satisfied in the form 
T = URq'^U^ R^ as (/3 7^ — in general. But the phase 
matching condition was derived assuming ajj <^ 1, and 
it cannot be used with \a) as the effective target state 
because {(t\U\0) = {0\WRtU\0) is close to 1. That is 
why we converted i?^ to R~'^. 

Now applying T on the initial state U\ffl rotates it 
towards the state |t) by an angle 26* sin ^ [7]. After n 
iterations of T, we get 

r"C/|0) = cos6l„|cr) +sin6'„|r), 6i„ ^ 9 [l + 2ns\u- 

^6) 

For n = [7r/(46'sinf )J, 6l„ is close to 7r/2 and T"f7|0) 
is close to \t). Further iterations of T rotate the state 
away from |t), displaying a cyclic motion in the two- 
dimensional subspace as in case of Grover's algorithm. 

To understand the significance of the state |r), we 
use the expansions U\0) = ^~Ujo\i) and /?|?7|0) = 



\j) in Eq. ([2]), and obtain 



sin 9 



cos tie 



Here cos9 = |(0|;7ti?t[7|0) | = J2 



smce 
1 



(7' 



we have the bound cos 9 > 



(7) 



and 



2al 



or 16*1 < 2au. Hence, (j|r) = UjoO{au) for 
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j ^ T, and the projection of |r) on the non-target sub- 
space is ^J^j^T 1 01''') P — 0{au)- This projection is 
very small, which makes |r) almost a state in the target 
subspace, |(i|T)| = 1 - O(a^) ^ 1. 

The number of queries needed to reach the state |t) is 
twice the number of iterations of T, as each iteration uses 
two queries. We therefore have Q = 7r/(20sin^). The 
normalization condition for Eq. ([7]), |(t|t)| = 1, gives 

sin^e^ (1 - cos 61)2 -t- ^4|;X,oP cos 6* sin^^. (8) 



For small 9, this yields 9 — y X^jsT ^l^jopsm^^, and 

Q = , (9) 

4sinfyE,,TlC^.oPsin2^ 

For later reference, we point out that the state |t) is 
close to target state only because 9 — 0{au)- That is 
true for the Rt transformations which act only within the 
target subspace, but may not be true for general selec- 
tive transformations 5*4 which perturb non-target states 
also. More generally, Eq. ^ provides |(i|r)| = O(^), 
and the iterative algorithm can amplify the projection on 
the target subspace by a maximum factor of 1/9. That 
may be too small for a general selective transformation 
to reach a target state. In section III, we use the idea of 
recursion to overcome this limitation. 



A. Comparison with Grover's algorithm 

When {S'f,S'o} — {/t,/o}i i-e. when ip ~ (f)j — tt, the 
operator T is simply 2 steps of the quantum amplitude 
amplification algorithm. To demonstrate the difference 
between T and for general {S'f , 6*0}, consider the sit- 
uation where St — Rf, i.e. rotation angle (j) is the same 
for all target states. Then T = UR^'^W R^'^'UR^W Rf 
while = UR'^U^RfUR^WRf. The quantum ampli- 
tude amplification algorithm succeeds in this case only 
if the phase-matching condition is satisfied, |0 — (/?! <C 
au [J] . On the other hand, there is no such restriction on 
our algorithm which succeeds using 7r/(4ac/ sin ^ sin ^) 
queries. As Grover's optimal algorithm takes TT/{Aau) 
queries, the slowdown is only by the constant factor 
l/(sin ^ sin -1). As long as (t>,ip are not very small, not 
much is lost, and hence almost any selective transforma- 
tion can be used for quantum search. 

This particular case has been experimentally verified 
on an NMR quantum information processor 8], which 
compares the performances of Grover's and our algo- 
rithm for small ajj and (j) ^ ip. The experimental data 
confirms the theoretical prediction that our algorithm 
succeeds in getting the target state while Grover's al- 
gorithm does not. In a more general case, the rota- 
tion angle 4>j can be different for different target states 



and St = Rt = X]j It is shown in an up- 

coming paper [§] that in this case, iterating the oper- 
ator Q — UR^U'^Rt amplifies only those target states 
which satisfy the phase-matching condition i.e. for which 
— <p| <C |tjo|- (If there are no such target states then 
iterating Q will not succeed in getting a target state.) 
The target state is obtained after 0{\/a'jj) iterations 

where a'jj = y^Ej:|0,-v|«|)7,o| l^^jol^ - ^u, and the al- 
gorithm suffers a slowdown by a factor au jc^jj- There is 
no such restriction on our algorithm and the full ampli- 
tude along the target states can be utiHsed, irrespective 
of any phase-matching condition. 

Quantum amplitude amplification is often described 
as a rotation in the two-dimensional space spanned by 
the initial state C/|0) and the target state \t). Here, we 
have provided a new insight suggesting that it is better 
to interpret quantum search as a rotation in the two- 
dimensional space spanned by the initial state t/|0) and 
the oracle- modified initial state i?|j7|0). T is then the 
fundamental unit of quantum search rather than Q . The 
advantage of the operator T is that it uses the selective 
transformations and their inverses in such a way that 
the phase-matching condition is effectively satisfied to 
produce a successful quantum search. 

B. New search Hamiltonian 

Grover's algorithm is a digital algorithm in the sense 
that it uses a discrete set of unitary operators and ap- 
plies them sequentially on the initial state to reach the 
target state. Farhi and Gutmann developed an ana- 
log version of the algorithm [l^], which shows that any 
initial state, when evolved under a particular search 
Hamiltonian for a certain amount of time, will evolve 
to the target state. Their search Hamiltonian is given 
by Ufg = Wt/|o) + W|t), where Hy|o) = / - [/|0)(0|C/t 
and 'H\t) — I — \t){t\ are projector Hamiltonians. More 
general search Hamiltonians have been presented subse- 
quently by Fenner [ll[ , and by Bae and Kwon . 

The algorithm developed above suggests a new search 
Hamiltonian 

'Wncw — Ti-mo) + W^jt^iQ^ — 'Hu\G) + R\'Hu\o)Rt- (10) 

The second term of Tincw is just the first term but in 
a basis rotated by the oracle transformation Rt- Tincw 
can be analysed the same way as was done by Farhi and 
Gutmann, in the two-dimensional subspace spanned by 
C/|0) andi?|c/|0). When evolved using Tincw for a certain 
amount of time, the initial state becomes the state |r), 
which is very close to a target state as shown. 

Tincw has certain physical implementation advantages 
over Ti-FG- Consider the situation when implementation 
errors perturb Hfg to (1 — s)H{j\Q'j + (1 + s)T-l\t), i-C one 
term is enhanced while the other gets reduced. Analysing 
this perturbed Hamiltonian as is done in Ref. [iSj, it 
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is easy to see that one reaches a target state only if 
|s| < 0{aif). This is analogous to the phase-matching 
condition, and as ajj ^ 1, it is a strict condition. There 
is no such restriction, however, on the new search Hamil- 
tonian as it is the sum of the same term in two different 
bases. For example, calibration errors remain effectively 
the same for both terms, making Tinow robust with s — 0. 



III. RECURSIVE ALGORITHM 

In this section, we consider those diagonal selective 
transformations {5t,S'o} which may also perturb the 
non-desired states, unlike the transformations {Rt,Ro\ 
discussed in previous section which leave them unper- 
turbed. We assume the perturbations to be small, i.e. 
\\St — It\\ = At and ||S'o — /oil = Ag where At and Aq are 
small. More explicitly. 



For such transformations, the iteration of operator 
USqW St on the initial state U\0) may not give us a tar- 
get state. As 0{l/au) iterations of UloWit on U\0) 
gives us a target state, it is easy to see that as long 
as {At,Ao} < 0{au), iterating the operator USoU^St 
on U\0) will also bring us close to target state. But 
when {At,Ao} > 0{au), the phase-matching condi- 
tion required for a successful quantum search may not 
be satisfied. Another way to see this is to analyse the 
eigenspectrum of USoWSt- Its distance from UloWit 
is 0(At, Ao). The two eigenvalues of UloWit, relevant 
for quantum search, are separated by 0{au). A per- 
turbation greater than 0{au) will in general shift them 
too much to maintain a successful quantum search. Note 
again that we are considering au <ti 1, which makes the 
iterative quantum amplitude amplification very sensitive 
to small errors. 

The recursive quantum search algorithm is defined, at 
the m*'' level, by the relation 



(12) 



with Uo = U. At the first level, ;7i|0) = {USoWSt)U\0) 
is a simple generalization of the quantum amplitude am- 
plification step GU\0). But at higher levels, the opera- 
tors Um involve {S't,S'o} as well as {5'|,S'o}, and can- 
not be expressed using repeated iterations of a single 
operator like G- (For instance, U2 — UiSoU\StUi = 
USoU'^StUSt)U^S\USlWStUSoWStU involves opera- 
tors 5* and 5*^ in a non-periodic pattern.) The idea of 
recursive quantum search is not new. It has been used 
by Hoyer et al. [l3| and by Grover [l^ for specific error 
models, as discussed in the next section. What is new 



here is the demonstration that recursion works even for 
general errors. 

The number of queries used at the m*'' level of re- 
cursion is determined by the relation ~ 3qm-i + 1, 
since qm~i is the number of queries used by Um-i and 
St needs one extra query. Using the fact that qo — (as 
implementing U does not need any query), we get 
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9(3"). 



(13) 



The recursive algorithm increases the number of queries 
in a geometric progression with the level number, a fac- 
tor of 3 in the present case. On the other hand, the 
iterative algorithm increases the number of queries in an 
arithmetic progression with the iteration number, a step 
of 2 in the algorithm of the previous section. We will see 
that the larger jumps in the allowed number of queries 
for the recursive algorithm are not a major disadvan- 
tage, because the total number of queries needed to ob- 
tain the target state remains about the same. (The worst 
case overhead is a tolerable factor of 3 in the number of 
queries.) 

At the first level, the initial state U\0) evolves to 
Ui\0), whose projection on the target subspace is au^ = 



X^jeT l(^i)io P- recursive quantum search, what 
matters is the amplification factor 



aui 
au 



E,eTl(t^: 



i)jo\ 



ST \ UjO 



(14) 



and the target state can be obtained using 0{l/a^^''^) 
queries. To get the nearly optimal algorithm, the ampli- 
fication factor AC should be as close to 3 (the number of 
U-type operators used by Ui) as possible. We will show 
that for small {A(, Aq}, k is indeed close to 3, and the 
performance of recursive algorithm is close to the optimal 
algorithm that takes 0{l/aij) queries. 

We estimate k by estimating the ratio pj = 
\{Ui)jo/Ujo\ for j e T. In terms of pj, we have 



(15) 



Clearly if pj is close to 3 for each j £ T, then k is also 
close to 3. To find pj, let 



|V)=5t[/|0)=^L/,oe^*^|j), 



(16) 



so that Ui\0) = USoW\ip). We decompose So as 5*0 

+E,^o e''^' |j)(j| leaves the |0) 



where5'^ = 

state unchanged but acts like "^0 on all the other states, 
and Rq" is a selective phase-rotation of the |0) state. We 
have Ui\0) = US'^WUR^^Wlij). With UR^^'Wltp) = 
\iP)-{l-e'^'>){Q\W\ip)U\0) and l-e'^" = 2e*''°/2 ^qs 
we get 



Ui\o)^us',u^m, 



(17) 
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where 



i/io/2 



COS ■ 



Mo 



P) \J)- (18) 



Here (3 = {0\W\^Ij) = \U,o\^e"f^ . As .^^ - tt^j^ + e, 
the bound \ej\ < At gives 

(1 - Re{f3)) < 0.5A2 + 2afj , \Im{(3)\ < A*. 



Since and afj are smah, we can write /? 
where |^| < At. Then 



(19) 
(20) 



where /3' = cos and = C - + f^- The bounds 
on P,£^, fiQ and give 

(1-/3') < 0.5A2 + 0.125A2 + 2a^, 

< 2At+0.5Ao. (21) 

Using Eq. we get \{jm/U,o\jeT = |1 + 2/3'e*«^|. 

The bounds on /?' and then yield 



01^') 



<lA't+lAtAo + lAl + 4alj. (22) 



Special Case: Consider the situation 5*0 = Rq , i.e. 
S'q = I. In this case, Ui\0) = \tp'), and we have pj = 
\{jW)/Ujo\ which obeys the bound ^ for j G T. Using 
Eq. ([Ill), we get 



(3-k) < ^A^ + ^AtAo- 



(23) 



Thus the projection on the target subspace is amphfied 
by a factor close to 3 as At, Aq and au are small quan- 
tities. The main idea behind recursion is to note that 
the above analysis holds for any unitary operator U, and 
hence it also holds for Ui which is a unitary operator. 
Therefore, U2 = UiSoUlStUi wiU obey 



(3 - K2) = ( 3 - ^ ) < -A^ + ^AtAo + ^Ajj + 4a 
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(24) 

where au^ — y^X^jer I0'l^2|0)p. Thus the projection 
on the target subspace is amplified again by a factor 
close to 3, making the total amplification close to 3^ = 9. 
Continuing the process, the ra*-^ level of recursion gives 
"(7™ = nlli i^io^u, where (3 - ki) < |A( + |AtAo + 



^Aq -|-4a^j_^. As long as 
fication factor obeys 3™ > YViLi — where 



As long as afj <C 1, the complete ampli- 
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A? - -AtAo 



1 



A2 
^0- 



(25) 



at least aur„ = 0(k'"q!j/). We can always choose m such 
that the condition afj^ = c ^ 1 is satisfied, and then 
repeat the algorithm times to get a target state. The 
number of queries required by the algorithm to get a 
target state is, therefore, at most Qm — 0{3^°^~^^^°'"^) = 

log— 3n 



0{l/ajj'^'' '^). In other words, the query complexity of the 
'log 3 



algorithm is ©(a^*-^"'^^-*), with 



<p = 



, _ - 1 < 0.71A2 + 0.20AtAo + O.IOA^. 

lOgK / 

(26) 

General Case: For more general 5o transformations, 
the state f/i|0) = US'^p^ip') is not equal to S'q is 

close to identity, however, and \\USqU^ — = HS'q — /|| = 
Aq. Upto a phase factor, we have 

{j\US'^U^ = cos7j(j| + sin-ij{xj\, 



(27) 



where \xj) is a normalized vector orthogonal to As 
\\US'f^W — /|| = Aq, we have the bound \ < Aq so that 
siwyj w 7j . Now 

(C/i),o - UlUS'oU^l^/) = cosjjiM') + lj{x,\ij'). (28) 
Using Eq. (|20p for \tp'), we find the ratio pj to be 



cos7je*'^^(l-|-2/3V«0 + 7j- 



(29) 



As IV') = UR^''W\tP), we have \{i^'\U\0)\ = \{ij\U\0)\ 
Hence, up to a phase factor. 



\^')^PU\0)+P\y) , P=^l^W, (30) 

where \y) is a normalised vector orthogonal to U\0). The 
bound on /3 (HH) implies P < y/Af + 4a\j. Eq. ([29]) then 
reduces to 



P]eT 




Csj 1 , 




Cij 


— cosjje'''^' (l - 






C2j 














(31) 



Since cosjj = 1 - O(A^) and 1 + 2f3'e'^'^ = 3 - 
0(A2,A2,AtAo) (as proved earlier), we have \Cij\ ~ 3 
for small {At, Aq}. Using the definition (|27p of {xj\ and 
the bound < Aq, 

{xj\U\0) ^ U,o- 



COS"fj 



i2\ 



= C/joO(Ao), 



(32) 



which makes — 0(Aq) and \Cij + C2j\ = 3 — 

0(At , Aq, AtAo). The ratio pj^r will then be close to 3 
iff 



This analysis shows that m levels of recursion can be 
used for amplifying the projection on target subspace to 



3jl 



< 3. 



(33) 
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By their definitions, the vectors \xj) and \y) depend upon 
the eigenvalues of Sq and St respectively. In most cases, 
the eigenvalues of these two different operators are un- 
correlated (in case they are correlated, we need to ran- 
domize one of them by random operations), and hence 
I a;) and \y) are two relatively random unit vectors in the 
A^-dimensional Hilbert space. So the expectation value 
of their inner product is and the above 

condition translates to 

^«3|f/,o|, (34) 

As long as this condition is satisfied for all j € T, the 
ratio PjeT and the amplification factor k are close to 3. 
More precisely, 

K = 3-0(A2,A2,AtAo). (35) 

If this condition is not satisfied for a particular target 
state J, then the amplitude along it will not be amplified 
by the recursive algorithm as if it were a non-target state. 

The condition ((M)) is only a sufficient, not necessary, 
condition for k to be close to 3. If it is satisfied for the 
first level of recursion then it is automatically satisfied 
for higher levels as \Ujo\ < \{Ui)jo\ for any I. Also, even 
if this condition is not satisfied then amplification may 
still be possible by a factor greater than 1, but not close 
to 3. Note that if or /? is 0(J7,o) then the condition 
is satisfied. It can be shown that this is the case when 
either of St or Sq becomes a selective phase-rotation Rt 
or Rq (the special case discussed earlier corresponds to 
Sq ~ Rq). Also, the condition is always satisfied for 
U &s Wjo ^ 1/Vn and jjP < Aq V'Af+I^ < 1. 

A. Comparison with Grover's algorithm 

When {St, So} = {/t,/o}, the recursive algorithm re- 
duces to the iterative Grover's algorithm and the optimal 
query complexity of 0{l/au) is achieved. The state at 
m*'* level of recursion Um\0) is nothing but (3™ — 1)/2 ap- 
plications of U laWlt on the initial state U\0). Explicitly, 
with ij = It and /q — Iq, 

u^+i = {uioU^it)'"-uioU\ilu4u^y"^it{uiountr'-u 
= {uiouUtY^-uioU\itUioU^Y"^it{uioU^ity"^u 
= [uhuntf^^+^u. 

With go = 0, = (C//of/Ut)(3'"-i)/2[/ is just quan- 
tum amplitude amplification, except for the jumps in the 
number of queries. 

In recursive quantum search, we are interested in the 
amplification factor n — au^jau of the projection on 
the target subspace, achieved by applying USoU^St to 
U\0). Detailed eigenspectrum of USoU^St is not of 
much relevance, since what matters is only one (rather 
than multiple) application of USqW St- In general, the 



state USoW'StUlO) = UIoWltU\0) + |A), where |A) 
has norm 0(At,Ao). k is certainly close to 3, when 
{At, Ao} <^ 0{au)- What we have shown above is that 
even when {At, Ao} 0{au), k can be close to 3. That 
is because what matters for n is not the norm of |A) but 
its projection on the target subspace, which can be small 
compared to ajj even when {At, Ao} 0{au)- 

The recursive algorithm needs 0{l/a^^^'^ ■*) queries, 
with A = 0(At,Ao) characterizing the size of errors. 
The increase in query complexity, due to nonzero A, is 
only a constant factor provided A = 0(-\/— 1/ logajy). 
This is a much better performance than the quantum 
amplitude amplification algorithm which needs A — 
0(1 /ajj) for success. Furthermore, the recursive algo- 
rithm can succeed even for larger A at the cost of more 
queries. 

IV. DISCUSSION 

Finally we consider the situation when {/g. It} are re- 
placed by non-diagonal operators {P,Q}. The iterative 
algorithm then evaluates {UP''U'fQ'<UPWQ)''U\0). Us- 
ing diagonal decompositions of {P, Q}, i.e. P = EpSoEp 
and Q = EqStEg with Sq and St diagonal, that becomes 

EgiVSlV^SlvSoV^StTVEllO), where V = E^UEp. 
The algorithm therefore converges to the target state in 
0{l/av) steps, provided {S't,<S'o} satisfy conditions for 
successful quantum search and (i?p)oo, {EQ)tt are close 
to 1. The condition (i?p)ooj {EQ)tt ~ 1 is important for 
any search algorithm, because only then we can rightfully 
call the transformations selective, performing an opera- 
tion on the intended state and leaving the other states 
alone. Thus, as long as T^o Uto, there is no significant 
slowdown in quantum search. 

Similarly, the recursive algorithm evaluates Um\0) = 
EQVmEl,\0) at the m*'* level, with 

K„ - Vra^lSoVi^^StVm^l. (36) 

As before, the algorithm succeeds, provided {St, So} 
satisfy conditions for successful quantum search and 
(£^p)oo, {EQ)tt are close to 1. 

Next wc point out a few applications of our algorithms. 

(1) Correction of Certain Systematic Errors: 
Quantum amplitude amplification is a repetitive applica- 
tion of the operator Q — UIqW It- Small errors in Q may 
accumulate over iterations to produce a large deviation 
at the end, causing the algorithm to fail. Completely 
random errors have to be protected against, using the 
techniques of quantum error correction and fault-tolerant 
quantum computation [l6j . That adds redundancy to the 
quantum states and gates, i.e. extra resources, to over- 
come small errors. For errors exhibiting specific struc- 
tures, however, it is worthwhile to investigate whether 
the dependence on quantum error-correction can be re- 
duced by designing quantum algorithms that are intrin- 
sically robust to these errors. 
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In this paper, we have studied a particular class of sys- 
tematic errors, those that are perfectly reproducible and 
reversible. For an imperfect apparatus in this category, 
we have presented two algorithms that exploit the struc- 
ture of errors and succeed in quantum search while the 
standard quantum search fails. These type of errors are 
not uncommon, e.g. the errors arising from imperfect 
pulse calibration and offset effect in NMR systems 8]. 
Thus our algorithms offer a significant flexibility in phys- 
ical implementation of quantum search. 

(2) Handling Errors in Workspace: The It trans- 
formation used in quantum search is implemented using 
an oracle. A typical implementation uses an ancilla qubit 
initialized to the ^"^^"'"^ state, and a C-NOT gate ap- 
plied to it from a Boolean function f{j). In general, f{j) 
has to be computed using the techniques of reversible 
computation, and has to be uncomputed afterwards to 
ensure reversibility. Inevitably, we need to couple our 
search-space to an ancilla workspace to implement /(, 
and the two get entangled. For a perfect algorithm, the 
workspace returns to its initial state at the end of the al- 
gorithm, and the search-space and the workspace get dis- 
entangled. But when there are errors, the workspace may 
not exactly return to its initial state, leaving some entan- 
glement between the search-space and the workspace at 
the end. That deteriorates the performance of quantum 
search, and our algorithms come to rescue in such cases. 

Let H = Hs Hw be the joint Hilbert space of 
the search-space and the workspace. The perfect or- 
acle is It = (|j)(j|/o>i <^ i-i) + lj)01/(j)=o ^ I)- 

In case of imperfect oracles, it may become Q — 
(|j>(il/(j)=i lj>01/0)=o «) B), where A,B are 

unitary operators. First consider the case B = I, i.e. 
the workspace remains unaltered for /(j) — 0. With 
the diagonal decomposition A = EqStEq, we have 

Q = Ej (|J>01/0>1®^Q^*4 + IJ')01/0)=0®^)- That 
is equivalent to Rt of section II, performing a selec- 
tive phase-rotation by (pk of the effective target state 
\j) f{j)=i<^\EQ{(j)k)) inH, where \EQ{(j)k)) is the eigenvec- 
tor of A with the eigenvalue e"^*-' . Our iterative algorithm 
would use f = UI^U^S\UI(p'^St, where U = U®I &uA 
Jq is the selective phase- inversion of |0) = |0) ® |0u,) with 
|0u,) the initial state of the workspace. As shown in sec- 
tion II, iterating T leads us to a state ® li/"), 
whose projection on the search-space is a target state. 
The number of queries depends on the eigenvalues of A, 
but it will be 0{l/au) as long as the eigenvalues are 
away from 1. The same result applies if the operator A 
is different for different target states. Note that this is 
a much relaxed criterion than the phase-matching con- 
dition which demands the eigenvalues of A to be within 
0{au) of -1. 

When B ^ I as, well, the iterative algorithm cannot 
take us to a target state and we have to use the recur- 
sive algorithm. The condition that \\St — It\\ should be 
small, restricts A to be close to — / (unlike the iterative 
algorithm, which allows a much wider range of A) and 



B to be close to the identity operator. For small er- 
rors in the workspace transformations, therefore, quan- 
tum search works and complete elimination of the entan- 
glement between the search-space and the workspace is 
not necessary. 

(3) Bounded Error Quantum Search: Our recur- 
sive search algorithm is similar to the quantum search 
algorithm on bounded error inputs by Hoyer et al. 
(labeled HMW henceforth), except that our error model 
is much more general. HMW considered computationally 
imperfect oracles, which provide the correct value of f{j) 
not with certainty but with a probability close to 1. For 
instance, if j is a target (non-target) state, the Boolean 
oracle may output 1 (0) with at least a probability 9/10. 
We have considered physically imperfect oracles, where 
the errors affect the unitary transformations correspond- 
ing to the oracle. In particular, the algorithm by HMW 
(see facts 1, 2 in section 3 of [l3|) uses fixed unitary trans- 
formations (5'o)hmw7 (<5'i)hmw (amplitude amplification) 
and E'hmw (error reduction), with (S'i)hmw replacing the 
oracle It- Our algorithm applies to the situation where 
these unitary transformations themselves contain errors. 
We have shown that as long as the errors are small, quan- 
tum search is possible. 

Indeed, the HMW error model can be reduced 
to our error model. The HMW oracle transforma- 
tion O computes the value of f{j) using workspace 
qubits and stores it in a qubit. It takes the ini- 
tial state ajb'>|0^->|0> to Ej ajlj)(VP7l'/'ji)|l) + 
■yr^^pj|i/'jo)|0)), where iV'jh), b € {0,1} denote the 
workspace states. The probability pj is at least 9/10 
if /(j) — 1 and at most 1/10 if /(j) = 0. Consider the 
operator Qo — OIqO^ SiO instead of only O, where Si 
inverts the states with last qubit |1) and Jq is the selective 
phase-inversion of the 1 0^)10) state. The operator Go ^s, 
an amplitude amplification operator, and its eigenvalues 
are e^^'^J with sin^O-j — pj [3]. Hence for /(j) = 1(0), 
the eigenvalues are close to —1(1). This is similar to the 
workspace error model discussed above, where \\St — It\\ 
is small. 

Moreover, if we assume that there are no errors in 
workspace transformations, our error model can also be 
reduced to the HMW error model. We simply attach 
a qubit to the workspace in the ^"^J^"'"^ state. A con- 
trolled St transformation takes the qubit to the state 
(|0) + e*'^^|l))/V2, where e^"^^ are eigenvalues of S't. A 
Hadamard gate then transforms the qubit to the state 
e»0./2(cos^|O) -ism^ll)). Since (/-^ =Trf{j) + ej with 
small \ej\, we obtain the HMW model. 

The difference arises when the workspace transforma- 
tions of the HMW model also suffer from errors. To get 
rid of these errors, we cannot keep on attaching extra an- 
cilla qubits till the new ancilla qubits are free of errors. 
Our results show that there is no need to worry about it, 
and recursion works as long as the errors are small. 

A peculiar feature of the HMW error model is that 
the imperfect oracle can be used to simulate an almost 
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perfect oracle by making 0(log N) oracle queries. There- 
after, the standard quantum search can be used. In our 
model, we cannot simulate It using St- In fact, we have 
shown that there is no need to simulate It] St is good 
enough for quantum search as long as it is close to It- 
More importantly, our algorithm also works when Iq is 
affected by errors, a case not considered by HMW. 

To conclude, we have presented two algorithms which 
allow a significant flexibility in the selective transforma- 
tions used by quantum search. The iterative algorithm 
takes 0{-\fN) queries and requires the oracle to be neu- 
tral for non-target states. But the oracle may mark the 
target states by phases other than phase-inversion, and 
hence almost any oracle transformation is good enough 
for quantum search. The recursive algorithm tackles 



the situations when the oracle perturbs non-target states 
also. For error size A, it reaches a target state using 
0{\/N ■ N'^'^'^ ^) queries. Needless to say, errors are 
inevitable in any physical implementation of quantum 
search. As long as the errors are small, the algorithms 
we have constructed are more robust and better adapted 
to physical implementation than the standard quantum 
search. 
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